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XVI. On the Times of Sudden Commencement of Magnetic 
Storms. By S. CHAPMAN. 


CoMMUNICATED BY Dr. C, Caren, F.R.S. 
RECEIVED Marcu 22, 1918, AND, IN REVISED FORM, AprRin 13, 1918. 


§ l. As at first communicated, this Paper was written in 
April 1917, or thereabouts.. Some revision has seemed de- 
sirable, partly in view of a subsequent change in my ideas 
concerning the nature of magnetic storms, and partly because 
a discussion of the same data by Dr. Chree,* with which I 
was not then acquainted, rendered certain introductory matter 
superfluous. 

In 1910 a Paper by Dr. Bauer} renewed interest in the 
question as to how far magnetic storms commence simul- 
taneously at different parts of the earth. W. G. Adams and 
W. Ellist had shown, long before, that the commencement- 
time is nearly the same, in great magnetic storms, over a large 
part of the earth ; as regards intervals of time amounting to 
afew minutes, however, the question was stillopen. Dr. Bauer 
claimed in his Paper to have demonstrated that the time of 
propagation “ round ” the earth was actually of the order of 
3 or 4 minutes, and he indicated the supposed direction and 
speed of propagation in the cases of several particular storms. 
These conclusions excited some controversy, being contested 
by various authors.§ In order to assist in settling the point, 
Dr. Bauer collected, with the aid of the directors of 32 con- 
tributing observatories, new data for the commencement-times 
of 15 magnetic storms.|| This set of material is the amplest 
and most suitable which has so far been available for the 
discussion of the subject, and there can be no doubt that its 
compilation has increased the interest and care taken in the 
accurate time-scaling of magnetograph registers. 


* C. Chree, ‘“‘ Proc.” Phys. Soc., XX VI., p. 137, 1914; also zbid., XXIII, 
p. 49, 1910; and “ Nature,” LXXXVL, p. 79, 1911. 

+ L. A. Bauer, “ Terr. Mag.,” XV., p. 9 (also pp. 107, 221), 1910; 
“ Nature,” LXXXV., p. 306, LXXXVL, p. 9, 1911. ; 

{t W. G. Adams, “ Brit. Ass. Rep.,” 1880, p. 201, 1881, p. 463; W. Ellis, 
“ Nature,’? XXIII, -p. 33, 1880, ‘‘ Proc.” Roy. Soc., LIL, p. 191, 1892 ; 
also van Bemmelen, ‘‘ Amst. Acad. Sci..”” Versl., XV., p. 250, 1906, Proe., 
1X.,_p. 266, 1906, Proc., 1908; also Birkeland, Norwegian Aur. Pol. Exp., 
1902-3, L, p. 63, 1908. 

§ Faris, ‘‘ Terr. Mag.,” XV., pp. 93, 213, 1910, and XVL, p. 109, 1911 ; 
Chree, I.c. ; Krogness, ‘‘ Nature,” LXXXV., p. 170, 1910; _ Walker, ibid., 
p. 236, 1910; van Dijk, ibid., LXXXVI., p.. 44, 1911; Birkeland, ibid., 
LXXXVL., p. 79, 1911. 

| “Terr. Mag.,” XVI., pp. 85, 163, 1911. 
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Up to this point the discussion seemed to have terminated 
without bringing about agreement between the opposing sides. 
The verdict of Dr. Chree and Dr. Krogness was that Dr. Bauer’s 
observational conclusions were not established, so that also the 
theory which he based upon them could not be maintained. 

No discussion of the subsequent collection of data has 
appeared in “ Terrestrial Magnetism,” and at the time of my 
own examination of it I was unaware that the material had 
already been discussed by Dr. Angenheister* and Dr. Chree.t 
So far as regards the above-mentioned. theory, the three dis- 
cussions agree in concluding that the new data afford no 
evidence in its favour, but rather to the contrary. In view, 
especially, of Dr. Chree’s very complete treatment it would 
be useless to labour this point further; as regards the dis- 
cussion of instrumental points, and the grouping of the time- 
data according to latitude and longitude, a general reference 
to his Paper will suffice. The justification for bringing forward 
this further discussion must rest upon the novelty of the point 
of view underlying it, which is different from that adopted by 
either Dr. Angenheister or Dr. Chree. Although no definite 
positive results are here derived from the existing data, the 
method of treatment does, I think, give promise of such when 
more material is forthcoming. 


§2. In the various Papers cited only three possibilities seem 
to have been contemplated regarding the commencement-times. 
Thus, Dr. Chree{ remarks: “As regards these ‘ sudden’ 
changes three things are conceivable : they may be absolutely 
simultaneous at different stations ; there may be a very small 
difference of time, corresponding to the rate of propagation of 
electromagnetic waves ; or, finally, there may be, as Dr. Bauer 
concludes, longer intervals, amounting to several minutes, for 
stations remote from one another. Under existing conditions 

‘of registration one cannot decide between the first two possi- 
bilities. As between these two and the third, a decision should 
not be impossible.” 

While reading Mr. Maunder’s Papers$ on the recurrences of 
magnetic storms at intervals equal to the synodic rotation- 


~ PEC “Nach. d. K. Ges. d. Wiss. zu Gottingen, Math.-Phys, 

+ Chree, l.c. (1914), 

t C. Chree, Proc. Phys. Soc., XXIII, p. 49, 1910. 

§ E. W. Maunder, ‘‘ M.N.R.A.S.,” Nov. 1904, April, May, 1905, Nov. 
1915; “Journ. Brit. Ast. Ass.,” XVI., p. 140, 1906; “ Astrophysical 
Journ.,” XXI., p. 101, 1905. 
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period of the sun (27-3 days), a fourth possibility naturally 
suggested itself. Mr. Maunder’s results suggest that storms are 
occasioned by some solar agent which is transmitted along 
narrow, well-defined streams, issuing from and rotating with 
the same period as the sun. The direction of solar rotation is 
such that the streams, if suitably situated, will overtake the 
earth on the afternoon or P.M. side.* Further, since the 
earth’s angular diameter, viewed from the sun, is only 18”-6, 
it is easy to calculate that the time necessary for a stream to 
sweep right across the globe is almost exactly 30 seconds. 

The above considerations suffice to explain how, without 
regard to any particular theory as to the precise mechanism 
of magnetic storms, the presumption arose that the relative 
time of commencement of a storm, at different stations, 
depends mainly on the orientation of the latter, at the time, 
relative to the sun. The determining factor, it is thus sug- 
gested, is neither mainly latitude, nor “ absolute ”’ longitude, 
but principally longitude relative to the sun. Now this relative 
longitude depends simply on the local timet at the station at 
the commencement of the storm. The local times of com- 
mencement, therefore, form the basis of classification of the 
data in the present discussion. 


§3. Again, the order of magnitude of the time-differences to 
be expected, between the commencements at different stations, 
is suggested by the time (30 seconds) taken for a stream to 
sweep across the earth. The differences, therefore, are to be 
reckoned in fractions of a minute, or in seconds, and not, as 
on the second and third of the hypotheses mentioned by 
Dr. Chree (§2) in small fractions of a second, or in whole 


* Tf we imagine ourselves looking from the sun towards the earth, the 
latter will be seen to travel along its orbit from right to left. The solar 
streamers rotate in the same direction, but, moving much more rapidly, 
overtake the earth. The earth and sun revolve in the same direction, so 
that the hemisphere visible from the sun (the ‘“ day ” hemisphere) will be 
turning from left to right. The local time for all stations to the right of 
the terrestrial meridian plane containing the sun will be after noon, so that 
the right-hand hemisphere will be termed the P.M hemisphere (post meri- 
diem). In this paper it will be assumed, for simplicity, that the axes of 
earth and sun are perpendicular to the ecliptic, or plane of the earth’s orbit. 
The earth will also be regarded as divided into four quadrants relative to 
the position of the sun at the time, the local time in the quadrants, which 
we may respectively term the night A.M., the day A.M., the day P.M., and 
the night P.M., ranging from midnight to 6 a.m., 6 a.m. to midday, midday 


to 6 p.m., and from 6 p.m. to midnight. 


+ Local time measures the longitude, at the rate of 15° per hour, from 
the ‘‘ midnight ” meridian, which is in the plane containing the sun, but on 


the dark side of the earth. 
Q 2 
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minutes. It is true that an electromagnetic change occurring 
at any station will produce some effect, propagated with the 
speed of light, at all other stations on the earth, but the 
magnitude of the effect diminishes so rapidly with increasing 
distance from the source that it soon becomes inappreciable. 
The initial impulse of a magnetic storm, like the diurnal 
magnetic variations, is mainly produced by atmospheric 
currents situated in the comparatively near neighbourhood 
of each station. The time differences of commencement at 
different stations are due to the lack of simultaneity of arrival 
of the solar agent of excitation at those places. 

Since this discussion was first completed, having carefully 
studied the nature of magnetic storms, my then vague ideas 
have given place to fairly definite conclusions regarding their 
origin and mechanism.* These do not enable me, however, to 
conclude any more definitely than at first whether the precise 
range of commencement-time to be expected, in the occurrence 
of a storm at stations distributed over the earth, is likely to 
be rather more or rather less than the above 30-second interval. 
The electric corpuscles which compose the streams do not 
impinge directly upon the atmosphere, but are deflected by 
the earth’s magnetic field, while still at distances comparable 
with the earth’s diameter. Prof. Stérmer ¢ has made calcula- 
tions as to the paths of single electric corpuscles in the earth’s 
field, and has found that the variety of possible types is very 
great. His calculations appear to be too ideally simplified, as 
yet, to give results accurately representative of the actual 
conditions, the mathematical difficulties even in the simple 
case of one electron being very formidable. The magnetic 
storms themselves indicate that the particles are precipitated 
more upon the P.M. than upon the A.M. hemisphere,} and as 
the former hemisphere is also the one on the side of approach 
of the streams, it is not unlikely that the initial impulse of a 
storm occurs first in that hemisphere. The discussion of the 
data of this Paper seems to favour this view, though the 
evidence is too slight to warrant a definite conclusion. It 
may be said, therefore, that at present theory gives but little 
guide as to the precise results to. be expected, and that more 


* These conclusions are described in a Paper communicated to the 
Royal Society cn April 10, 1918. 

+ Stormer, ‘‘ Arch. d. se. phys. et nat. Genéve,”’ 1907, and various articles 
in ‘“ Terrestrial Magnetism ” and elsewhere. 

t Cf. the footnote-to p. 207; also the Paper on magnetic storms re- 
ferred to in the last footnote but one. 
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observational material is required in order to ascertain the 
actual facts. 


§4. The data to be considered are confined to the times 
of commencement as measured from the horizontal force 
record alone. The sudden perturbation at the beginning of 
a storm is generally most pronounced, and therefore most 
accurately measurable, in this element. The principal source 
of error in the data appears to be the uncertainty as to which 
is the precise movement on the registers which marks the 
sudden commencement. If proper precautions are taken, the 
time of a given movement on the record is measurable with an 
accuracy of about a minute. An inspection of the data 
suggested that (with perhaps one exception) the systematic 
error of time determination, at any of the 32 observatories 
which contributed data, did not exceed one minute. I there- 
fore decided to reject any observation which stood out from 
the mean time of commencement for the particular storm by 
two minutes or more, on the ground that the difference pro- 
bably arose from a choice of the wrong movement for 
measurement. Some observatories gave more than one com- 
mencement-time for one or two of the storms, accompanied 
by drawings showing the magnetic movements to which the 
times corresponded. In these cases it was usually possible 
to choose, with little uncertainty, the precise movement 
which had been taken by the majority of the other observa- 
tories. The total number of rejected values was not large,* 
while the number retained was 342 ; if all the 32 observatories 
had given data for each storm, the total number of values 
would have been 480, but some observatories did not furnish 
complete data. The table (facing p. 212) shows the rejected 
values as well as those retained. 

The mean time of commencement for each storm was 
determined from the values given by the several observatories, 
excluding those which were rejected. The residual differences, 
“ actual value—mean,” are given in the table, for all the con- 
tributing observatories. Where a query is attached to an 
observation, this was done by the observatory which supplied 
it, and these values have in all cases been bracketted. A few 
other particulars regarding the table are indicated in the 
notes appended to it. The observatories are arranged in order ’ 


* Counting those rejected on ground of magnitude alone, there were 
31 negative and 22 positive residuals bracketted. 
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of longitude, and the storms in order of date. Columns are 
added giving the systematic difference of time from the mean, 
and also the mean numerical residual, for each observatory ; 
in a row beneath, the mean numerical residual for each storm 
is given. The number of residuals concerned in each mean is 
indicated by a suffix. The unit in which all these* quantities 
are expressed is one-tenth of a minute. 

For some storms the data are clearly much better than for 
others ; this is probably owing to the greater definiteness of 
the initial perturbation in the former cases. Magnetic storms 
show marked differences in this particular. Storm 1 is an 
example of a “ good” storm—only two values are rejected 
out of 27, and the mean numerical residual of the other 25 is 
0-4 minute. The numerical mean residual derived from the 
342 values all taken together is about 0-7 or 0-8 minute. 
This should allow a clear indication to be given of times of 
the order of two or three minutes, such as were suggested 
by Dr. Bauer ; this, it will be seen, is not shown, the re- 
siduals, which are of about one minute in magnitude, largely 
consisting of accidental error. When it is remembered that, 
with the ordinary time-scale, 1 mm. on a magnetograph sheet 
represents about three minutes of time, this degree of acci- 
dental errcr is not surprising. 


$5. In accordance with the ideas described in §§ 2, 3, the resi- 
duals for each storm were divided up into four groups, accord- 
ing to which of the four quadrants mentioned in the footnote 
to § 2 the corresponding station was situated in at the time of 
commencement. This was done by considering the longitudes, 
west of Greenwich, of the four principal meridians with respect. 
to the sun, viz., the midnight, midday, 6 a.m. and 6 p.m. 
(local time) meridians. Thus, in storm 1, the commencement 
-oceurred at 19h 56-7m G.M.T., so that the midnight meridian 
was this amount (which in are is 299°-2) to the west of 
Greenwich. The longitudes of the four meridians (0, 6, 12 
and 18 hours local time) were therefore 299°-2, 209°-2, 119°-2, 
and 29°-2. In the table dividing lines representing these 
longitudes were drawn for the first storm (and in a similar 
fashion for the other storms), these lines being respectively 
light and continuous, broken, thick and continuous, and 
broken. The stations in the night A.M. quadrant at the 


* At the foot of the table, however, are other quantities, described in §5, 
which are expressed in units of 1 second. : 
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commencement of any storm are consequently those above’ 
the light continuous line and below the consecutive (broken) 
line. In this way it becomes an easy matter to group the 
residuals in the table according to their corresponding quadrant 
or hemisphere. 

At the foot of the table rows are to be found which give for 
each storm the means of the residuals in the four hemispheres, 
day, night, P.M. and A.M., and also the differences between 
the results for complementary hemispheres. These values are 
all given in units of 1 second, not 0-1 minute as in the re- 
mainder of the table. The mean difference for either pair 
of hemispheres in the case of no storm exceeds 60 seconds, and 
these mean differences for individual storms are clearly still 
affected by considerable accidental error. Any systematic 
quantities underlying them must, therefore, be of the order 
of 20 seconds or less. Taking all the residuals in any one 
quadrant together, we obtain the following results. The 
number of residuals, their sum* (in minutes) and their mean 
(in seconds), are given in order :— 


Night P.M. quadrant (84): Sum —5-0m., Mean —3-6s. 
Day P.M. - (89): ,, —3:8m., ,, —2:6s. 
Night A.M. 5 (SH). —l-Im...<,. —0 8s. 
Day A.M. (82): ,, -+7:2m., ,, -+5:3s. 


The quadrants are arranged in the numerically increasing 
order of their mean residuals. This order is not to be regarded 
as the true order of arrival of the initial impulse of a magnetic 
storm at the different quadrants ; more data would be needed 
to establish any such conclusion. There is perhaps enough 
evidence to favour the opinion that the P.M. hemisphere is 
the one first affected by the storm, though even this is doubtful. 
If we group the quadrants into pairs of hemispheres, we obtain 
the following results :— 


P.M. residuals (173): Sum —8-:8m.: Mean —3-2s. Snes 
agg PF Gilman ee ft 


Day residuals (171): Sum +3-4m.: Mean +1-2s. vend aces 
Night Met. on | 2a pe ee 


The latter result is certainly more unlikely than the former, 
Dut cannot be dismissed as (theoretically) incredible. In order 
to test how far accidental error might account for these small 

* It may be noted that the sum of all the residuals is —2-7m., which arises 


from an accumulation of small errors, less than 0-1m.,in the adopted mean 
times of commencement of the storms. 
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differences of time, I made the following two quite arbitrary 
subdivisions of the data :— 


Difference 
0-4s. 


Sum —1-:9m.: Mean —0-7s. 
Second 16 observations (170 residuals) : 


First 16 observatories (172 residuals) : | 
Sum —0-8m.: Mean —0-3s. 


’ 


Upper right-hand diagonal half (177 residuals) ‘) : 
Sum +2-5m.: Mean +0°8s. Difference 
Lower left-hand diagonal half (165 residuals) : 2-78. 
Sum —5-2m.: Mean —1-9s. 


The latter, which is about equal to the difference between 
the day and night residuals, confirms the unlikelihood of the 
latter being real. The matter for surprise in all cases, however, 
is the extreme smallness of the residuals, and this suggests that 
it should not be very difficult, even with present instrumental 
means, to determine any systematic time differences between 
the different quadrants, if these amount to so much as 10 
seconds. The extreme range (30 seconds, more or less) would 
not, of course, show in the mean of whole quadrants or 
hemispheres ; and the fact that the majority of observatories 
contributing data are in high latitudes likewise tends to 
minimise the differences to be expected. 


. ABSTRACT. 

The Paper is a discussion from a new standpoint of the data, col- 
lected by Dr. Bauer, for 15 magnetic storms. Maunder’s work on 
the recurrence of magnetic storms at intervals equal to the rotation 
period of the sun suggests that storms are due to some solar agent 
transmitted along narrow well-defined streams, issuing from and 
rotating with the sun. This suggests the view that the relative time 
of commencement of a storm at different stations depends mainly 
on the orientation of the latter at the time, relative to the sun, i.e, 
on the local time at the station, This forms the basis of the classi- 
fication in the Paper. ; 


DISCUSSION. 


Dr. Curer said that after hearing the Paper he was uncertain whether 
Dr. Chapman did or did not believe that the figures proved anything. 
He (Dr. Chree) did not think that they did. To deal with a question of a- 
few seconds’ difference at different stations, it would, he thought, be 
desirable to employ not merely a more open time scale, but magneto- 
graphs of greater and desirably -uniform sensitiveness. When a magnetic 
change became recognisable depended on the size of the movement and 
the sensitiveness of the magnetograph. An apparent difference in time 
between day and night hemispheres would naturally arise if movements 
tended to be larger in the one than the other. One objection which he 
had urged against Dr. Bauer’s views also applied to the hypothesis con- 
sidered by Dr. Chapman, viz., that as soon as a disturbance began at any 
part of the earth’s atmosphere it would naturally be propagated b> 


| a a a a ee ee 


Index number of storm, wit 
West , i| 2 3 Br mets F 
Observatory. longi- eo ae ARREST) ONG Haeck - ee ay, 
tude. . 1906, .1906, 1906, 1907, 1907, | 1907, 1908, 
July 29. Aug. 7. | Deo. 21.| Feb. 9. July 10.) Oct. 13. | Mar. 26. | 
19h56-7m 13h38 1m)21h30-6™) 14b13-0m)]4b99-7m| 7h43-9m 17b4] +4) 
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Falmouth* ...... 51 50-1,, 3 —l (34) (20) 13 18 
ce — | ete || 
Pilammias eere ence. 63-9 | 31-78 see ae —5 (25) 9 =6 
Porto Rico...... 65-4 18-1 N =2 —14 (34) 1 =s) 2 
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Agincourt ...... 79-4 | 43-7 ,, —4 (—20) 4 —5 —2 0 
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Batavia ......... 253-2 | 62s —6 | —18 —6 | —10 | —4 | (24) 
PoOungo00 wes coe. 263-5 18:9 nN 6 —3 Si Us 18 eae | 
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Wecledntieetesees 3556 50-8 ,, 3 14 (23) ad: = (36) | 
Val Joyeux...... | 358-0 48-8 ,, —7 —12 —6 —16 3 (—46) 
Biro. ' 359-5 | 40-8,, (39) =o \y(—28) | 
Mean numerical residual ............ 45. 853 650 Tao 656 Bea 
a ———_—__—. : - . — 7 = 
+ - ( Mean P.M. residual .| —7 45 3h St ra 5 
ie ame of | Mean A.M.residual.| 6, | —10, | —5, | —8. | —8,' i 
SiSoVablocn aM AIMS miececccsoes 13 —l4 —8 —4 —10 —4 
l second (not ) «, rier ‘ bie = 
0-1 seo are “ Day residuals | —10, 24,4 13, —2711 ho Si4 
as above) Night ” residuals .| —1,, | —35, —545 21, —4,, —9, 
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* These o bservatories gave times in whble minutes only. For ¢ 
t These times were in error by positive amounts varying from [ 
t A small movement, of perhaps local character, nearly coincid 
§ Several observatories gave data for two movements in this « 
one is here considered, which necessitates exclusion of dé 
To face p, 212, 


mean time of commencement. 
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e decimal 0 has been supposed added. 
8. 
scured the precise commencement in this case. 

t+ one being at 8h. 30m. approximately. Only the later 
ratories which measured only the former one. 
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electromagnetic waves to other parts. He was not aware of any direct 
evidence confirmatory of Dr. Chapman’s statement that the atmospheric 
currents causing the diurnal variation were mainly situated in the com- 
paratively near neighbourhood of the station. Also, even if this were 
the case, he did not see that any inference could be drawn as to what 
happened in the case of so different a phenomenon as a “ sudden com- 
mencement.’’ We did not even know whether the electrical currents 
causing the two phenomena were at the same level in the atmosphere. 
“ Sudden commencements ”’ were sometimes large, 50 or even 100¥. 
They were not instantaneous changes, the rise normally seen in hori- 
zontal force in ordinary latitudes taking usually four or five minutes to 
attain its full value, a curious feature being that the high value was 
generally retained for some time, sometimes, in fact, for several hours, 
especially at stations in certain latitudes. A consideration of the exact 
nature of the phenomenon was a desirable prelude to any theorising. On 
the historical question which Dr. Chapman had raised, as to the discovery 
of a tendency to “ repetition ”’ in magnetic storms, it was very natural 
for him as a Greenwich man to emphasise the really valuable work which 
Mr. Maunder had done in this connection, but it was only fair to remem- 
ber the much earlicr work of Broun, and the independent work of A. 
Harvey and the late Prof. Birkeland. One reason why a belief in the 
“repetition ” of magnetic storms had advanced so slowly was that sun- 
spot theories had been so often suggested by cranks that they were 
naturally somewhat suspect. In the case of magnetic disturbances, and 
Mr. Maunder’s Greenwich lists were no exception to the rule, the great 
majority did not have a ‘ sudden commencement,” and opinions as to 
when the average storm began or ended might differ by several hours. 
It was thus in general impossible to assign an exact value to the interval 
between two storms. Further, what one man called a “ storm” another 
did not, so where one might see a “' repetition ” another would not. A 
further complication was the existence of a diurnal period in disturbance 
fata, representing partly statistical imperfections and partly a real 
natural phenomenon. Dr. Chapman did not seem to realise that the 
27-day period was iust as manifest in quiet as in disturbed conditions. 
It was also important-to remember that it had presented itself in years 
like 1913, when sunspots were almost non existent, just as decisively asin 
years of many sunspots. 

Prof. Nrwaru said that when Maunder’s results were first brought for- 
ward it was felt to be a difficulty that they should involve an almost 
geometrical recurrence ; because, of course, the disturbances on the sur- 
face of the sun were travelling with varying periods. When the slide 
showing the periodic recurrence of certain storms was first shown by 
Maunder it was noticed that a 1 per cent. difference in the adopted 
rotation period would change the direction of the lines which indicated 
eoincidence of period to 45 deg. from the vertical. One saw on looking 
at the slide that these were many cases in which lines at different angles 
would suit. No one familiar with eclipse phenomena could doubt, 
however, that something originated in the sun. Many of the streamers 
sent out from the corona seemed quite straight, and it was difficult to 
resist the idea that sometimes the earth would be at the other end. Had 
Dr. Chapman considered the difficulty raised by Kelvin of the necessary 
energy relations ? The energy of the magnetic storm itself might be 
attributed to something abstracted from the energy of rotation of the 
earth by the particles after entering the atmosphere. 

Father Corrie thought Maunder was the first to connect the 27-day 
period with the solar streamers. There seemed to him to be a diffienlty 
about these streamers of charged particles. Schuster showed that 
owing to the mutual repulsion of the charges the stream would be dis- 
persed. He did not think that any particular storm was associated with 
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a particular sunspot, but with a disturbed area of the sun which might 
extend many degrees. His idea was that the dispersed streamers from 
the disturbed regions formed clouds of particles and that the earth may 
enter such a cloud, giving rise to magnetic storms. From the table it 
seemed to him that there was a slight preponderance of effect on the night. 
side relative to Greenwich—i.r., over the Pacific ocean. Would the 
presence of a large expanse of water be likely to affect the phenomena ? 

Mr. T. SmirH asked why the times of sudden commencement were 
chosen. He would imagine that similar time relationships would hold 
for any prominent feature. One would normally expect greater irregu- 
larities just at the beginning of the storm than when it was well under 
way. When it was known that a storm had commenced, wide scale 
runs could be started and measurements made on some subsequent out- 
standing feature. 

A Mremeer asked what irregularities might be due to lag in the different. 
instruments used in the different observations. 

Prof. Lrrs said that when Dr. Chree first brought the matter before 
the Society he had set some of his students to work out the correlation 
factor of Bauer’s figures. The factor was so low that there was evi- 
dently no basis whatever for Bauer’s theory. 

Dr. CHAPMAN, in reply, said he did not put any reliance on the actual 
figures given in the Paper. The observational error was too great. It 
was the method of treating the observations that he thought was of im- 
portance. He thought diurnal variations were produced at lower 
levels of the atmosphere than storms. By “‘ near neighbourhood ”’ he 
meant within about 1,000 kilometres or so. He did not think Maunder 
connected the storms with actual sunspots, but rather with disturbed 
regions. Of course, any proper motion of the disturbed region on the 
sun would alter the period in particular cases; but in examining a long- 
series it was best to take the synodic period since the others would be 
distributed on either side of this. He had not yet gone into the energy 
question, but did not think that this would present any insuperable 
difficulties. He did not think the streams would diffuse very much. He 
thought they consisted of a particles. They did not necessarily proceed 
radially from the sun, but might emerge in all directions from radial to 
tangential. This, however, would not seriously affect the time taken 
to pass across the earth. There was a definite reason for chosing the 
commencement times of the storms. When the storm is under way the 
magnetic state of the earth is fluctuating, and it is difficult to recognise- 
accurately particular features. There was certainly a lag in the measure- 
ments which differed in different instruments ; but the effect of this was 
greatly reduced in the present method of grouping the observations, since 
the same instrument on different occasions contributes results to different. 
groups. In this respect the method has a great advantage over classi- 
fication on a geographical basis, in which case any instrument is always. 
in the same group. 
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XVU. The Entropy of a Metal. By H. Stanuey ALLEN, M.A., 
D.Se., University of London, King’s College. 
ReEcEIvEeD Apnit 127H, 1918. 

§ 1. According to the theory of Ratnowsky* the entropy 
of a gram atom of a substance in the solid state may be ex- 
pressed in the form, 

( u £8 ; 
s= 4] Sede —x log (l—e I, 
0 


x3 1 


where N is Avogadro’s constant, k is the gas constant for 
a single molecule, and 2=hv,/kT=fv/T=0/T, h being 
Planck’s constant, vm the maximum vibration frequency of 
Debye’s theory, and © a temperature characteristic of the 
substance considered. 

Tn a previous communication f I have shown that the correct 
form of the approximation for smal! values of x, corresponding 
to high values of the absolute temperature 7, is 


Ay {4 x 
AE log t+ 46 bors 


where 3Nk=C,.=5-96 calories. 


§2. Up to the present time no direct evidence as to the vali- 
dity of this equation has been brought forward, but an im- 
portant Paper by Lewis and Gibson ¢ on the entropy of the 
elements now provides the data required for testing the 
formula. The work of these authors is based on the heat 
theorem of Nernst—the so-called third law of thermodynamics 
—which may be stated in the form: “In an isothermal 
process involving pure solids and liquids, the change in entropy 
approaches zero as the temperature approaches the absolute 
zero.” Following Planck, the principle may be re-stated 
in the more general form: ‘‘ The entropy of every actual 
substance in the pure state is zero at the absolute zero of 


temperature.’ Hence 
7dQ_ 'CaT 
a ad a 


* Ratnowsky, “ Deutsch. Physikal. Gesell. Verb.” Vel. XVI., p. 1033° 
4 


1914. 
+H. S. Allen, ‘‘ Proc. Phys. Soc.” Vol. XXVIII., p. 302, 1916. 


+ Lewis and Gibson, ‘‘ Am. Chem. Soe. Journ.” Vol. XXXIX., p. 2554, 
1917. 
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where C is the atomic heat. The absolute value of the entropy 
may be calculated, either at constant volume or at constant 
pressure, for any assigned temperature, provided C is known 
as a function of the temperature. The primary object of 
Lewis and ‘Gibson was to determine the entropies of the 
elements in their standard states, at the standard temperature 
of free-energy measurements, namely 25°U. or 298°K. For 
many substances the atomic heat at constant volume is given 
by the equation 
C.=} (L / 9), 


where / is the same function for different substances, and 6 a 
characteristic constant for each substance. Instead of making 
any assumption as to the exact form of the heat-capacity 
equation, the authors employed a graphic method of calculating 
the entropy of solid substances for which the equation has the 
same form. It is to be noted that 6 is defined by them as the 
temperature at which C, is one-half of the Dulong and Petit 
constant, so that when 7=6, C,=3R/2=} C,. 


§3. The following table contains the values of the character 
istic constants and the entropies for seven metals which have 
been made the subject of accurate investigations at low tem 
peratures. The column headed @ gives the characteristic 
temperature as defined by Lewis and Gibson, 9 is the character- 
istic temperature employed by Debye.* The entropy in 
calories per degree at 298° K., under the condition of constant 
volume, is given in the last two columns of the table. The 
observed values are quoted from the Paper by Lewis and 
Gibson,t whilst the calculated values have been obtained 
from the approximate form of the formula of Ratnowsky. 
The agreement between the observed and the calculated 


* It is noteworthy that the values of © in the table are almost exactly 
four times the corresponding values of §. It is interesting to test this 
relationship assuming Debye’s expression for the specific heat, which may 


be written 

Cy, 12 [" Bde 3x 

C, «J oef—1 e™—1 
I have to thank Prof. J. B. Dale for valeulating the value of the right hand 
side of this equation when x=4, that is when O=47. It is found that 
Cy/C, =0-503059, which differs from one-half by less than one per cent. 
This result shows that when Debye’s theory is employed, the characteristic 
temperature, 9, is almost, but not exactly, equal to four times the tem- 
perature at which Cy is one-half of the limiting value Co. 


+ These authors also give values for the entropy of the elements under 
constant pressure. 
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values is remarkably good. The letter a after an observed 
value indicates that the value in question is hardly likely to 
be in error by more than one-third of an entropy unit. 


ENTROPY ( OF THE Meats AT 298° K 


Plemcnt, 0. | 6. Se see re calculated. 
Aluminium ...... 95-5 5 3825 6:73 a 6-72 
WWOVS va citar 99-1 385 G 6-54 6-68 
Copper 78:2 315 K 791 a 7-79 
ANIC eae eencesteosen c 57-6 230 G, N 9-60 a 9-59 
Silyersiccey cosas ss 53:7 215 G 10:00 a | 9-98 
Cadmium bed “4 168 G 11-38 11-42 
WGA CH cc cectoesius 2-0 88 K, 8 15-lle 14-89 


G= Griffiths, K=Keesom and Onnes, N=Nernst, S=Schwers. 


The thermal behaviour of sodium * is somewhat exceptional, 
and for that reason it has not been included in the table. 
At low temperatures Messrs. Griffiths find 0 = 180, which 
yields as the calculated value of the entropy 12-21 calories 
per degree, whilst the value deduced by Lewis and Gibson 
from Dewar’s measurements is 11-43. 

Another case in which it is possible to institute a comparison 
between the theoretical and the experimental value for the 
entropy is that of solid mercury at the temperature of the 
melting point. Nernst and Lindemann found 0=96-6, so 
that at the melting point (7=234-1° K.) the entropy is 13-26 
calories per degree. The corresponding value deduced by 
Lewis and Gibson from the observations on the specific heat 
is 13-31 calories per degree. Here the agreement could 
scarcely be improved upon. 


$4. When z is not too large, the value of the entropy can 
be calculated from the formula of Ratnowsky expanded in 


the form 
xt xo ) 


S=Co = log “ag ‘3240 108864 |” 


The series is convergent when z is (numeraiclly) less than 27, 
but converges slowly when a is as large as 4, For large values 
of x the method of expansion given by Debye f may be em- 
ployed, or by combining the expression of Ratnowsky for 


* Bh. H. Griffiths and E. Griffiths, “ Phil. Trans,” A. Vol. CCXIV., 
p. 319,1914. Eastman and Rodebush, vu Journ. Am. Chem. Soc.”’ Vol. XL., 


p. 489, 1918. 
+ Debye, ‘‘ Annalen der Physik.” Vol. XXXIX., p. 789, 1912. 
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the entropy with that of Debye for the specific heat the entropy 
may be found from the formula, 


x 


S=40,+-C* E= —log (1 -e*) |, 


in which the values for C, calculated by Debye may be 
inserted. 

In the diagram (Fig. 1) the value of the entropy in calories 
per degree, calculated as described above, has been plotted 
against the value of 1/x=7'/0. For any assigned value of 0 
this gives an entropy temperature diagram. In the same 
figure are given the points (marked x) determined from the 
values of S, recorded by Lewis and Gibson as corresponding 
to selected values of 6/7. In plotting these points it has been 


“0 0:2 04 O06 08 10) Ike, 14 16 18 2:0 
Vx=T/O . 
Fic. 1.—VartaTIon oF ENtTRopyY WitH TEMPERATURE. 

©) Formula of Ratnowsky. <x Lewis and Gibson. 


assumed that O=4-0210, which is the relation given by 
Debye as corresponding to C,=3C.,. It will be seen that 
the two sets of points lie on a smooth curve, showing the close 
agreement between the theoretical and the experimental 
results throughout the whole range of tempe-atures con- 
sidered. 


_§5. The above results prove that the formula of Ratnowsky 
gives values for the entropy of a solid in very close agreement 
with those obtained directly from observations on the specific 
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heat at various temperatures. It therefore becomes a matter 
of special interest to examine the principal assumptions made 
in the deduction of the formula. These appear to be four in 
number. 


1. The whole internal energy, U, of the gram atom is made 
up of two parts, the “ vibrational energy,” #, and the potential 
energy /(v), which is a function of the volume. Even at the 
absolute zero of temperature the substance possesses a certain 
amount of potential energy. 


2. Beltzmann’s equation for the entropy is assumed, so 
that 
S=k log W+C, 


where W is the number of “ complexions ” of the system when 
it has a given energy. The constant of integration, C, is 
ignored or included as a constant multiplier in the quantity 
W, so that the equation reduces to 


S=k log W. 


This, as Planck points out (“Heat Radiation,” § 120), is 
equivalent to assigning a definite absolute value to the entropy 
S, and leads necessarily to Nernst’s heat theorem. 


3. It is assumed that 
eck 
oH T’ 
This implies a differentiation at constant volume. 


4. Following Debye, it is assumed that the number of 
vibrations within a frequency interval, dv is given by 


9Nv7dv/v3,,, 


where v,, denotes the maximum vibration frequency. It is 
known that Debye’s theory must be regarded as giving only 
a first approximation, since the existence of a sharply defined 
maximum frequency is a somewhat arbitrary assumption. 
Nevertheless the formula deduced for the atomic heat gives 
a very close approximation to the experimental results. 


§6. We may conclude that the hypotheses assumed by 
Ratnowsky in dealing with the entropy of a solid are probably 
justified as being at least approximately true. It is necessary 
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to point out, however, that the results of the present inquiry 
have no bearing on the further assumption made by Ratnow- 
sky that his expression for the entropy may be applied to the 
liquid state. This is a question which I have discussed in 
a former Paper dealing with the latent heat of fusion of a 
metal (loc. cit.), but it cannot be definitely decided till further 
information is available as to the applicability of Debye’s 
theory to liquids. 


ABSTRACT. 


An expression for the entropy of one gram atom of a substance 
in the solid state has been given by Ratnowsky. In a communica- 
tion to the Physical Society in 1916 the author gave the correct 
form* of the approximation required for high values of the absolute 
temperature in terms of Bernoulli’s numbers. The data required 
for testing the formula have been supplied in a recent Paper by 
Lewis and Gibson, who have given values for the entropy of the 
elements under the condition of constant volume, and also under 
constant pressure. These values were deduced from observations. 
on the specific heat assuming the truth of the heat theorem of 
Nernst, that the entropy of every actual substance in the pure state 
is zero at the absolute zero of temperature. It is found that the 
formula of Ratnowsky gives values for the entropy of a solid in very 
close agreement with those obtained by Lewis and Gibson. The. 
hypotheses assumed in the theory of Ratnoswky are discussed, and 
the conclusion is drawn that these are probably justified as being at. 
least* approximately true. 
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XVIII. On Tracing Rays Through an Optical System. (Second 
Paper.) ByT.Smiru, B.A. (From the National Physical 
Laboratory.) . 


ReEchivep Aprit 22, 1918, 


In the earlier communication* bearing this title the author 
gave formule by which a ray in three dimensions could be 
traced through a system of coaxial spherical refracting sur- 
faces. These formule are algebraic, and do not require the 
use of any tables ; in this respect they are in marked contrast 
to those generally used. The new formule are primarily 
intended for use in conjunction with a calculating machine ; 
the older trigonometrical formule have been arranged for 
logarithmic work. Experience has shown that logarithmic 
calculations take several times as long as the corresponding 
mechanical operations. 

The present Paper contains the modification of these 
formule generally used for rays in an axial plane, and, for both 
two and three dimensional cases, new formule which are 
universally applicable, so that the calculations could, if 
required, be carried out entirely by mechanical means. In the 
three dimensional case a modification is only required for the 
preliminary calculation described in the previous Paper. The 
present Paper also-includes the method followed in calculations 
relating to transverse focal lines. 

By far the greater number of rays that have to be traced 
through systems of lenses lie in a plane containing the axis. 
The system used in the previous Paper can be modified for 
these cases. Let 9, 9’, be the angles of incidence and 
emergence at the surface which separates media of refractive 
indices f,—-1 and yy. The inclination to the axis of the ray 
between surfaces Zand J+1 is y,; 4 and af are the distances 
between the vertices and the centres of curvature respectively 
of these surfaces. The radius of curvature of the surface 
Zis 7, and its curvature is Ry\=1/r,. The length of the per- 
pendicular from the centre of curvature of this -urface on to 
the incident ray is h,f, and on to the refracted ray is hy. To 
find the position of the final emergent ray the formule 


* Proc. Phys. Soc., Vol. XX VII., page 502. , d 
+ There are minor alterations from the notation employed in the previous 


Paper. 
VOL. XXX. R 
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Ay=h’) -1+@)_1 sin Wr-1 Pr sie ece> iar CA x (1) 
PP aller any A ee 
Bin Q. =A Rj. . ee RE at Sen) 
Bin op ys=/yhy 3. ys 29 ee 
Pr=Qn—tprqs (5) 
A, =b —I+N=1 (6) 


are successively applied at each refracting surface. 

Of these formule, (1), (2), (3), (4), (6) are equivalent to some 
of the general formule contained in the previous Paper. 
Equation (5) constitutes a modification which applies only to 
the particular case of rays in two dimensions. The formule 
as given above involve reference to tables of natural sines, 
since equation (5) refers to angles, and not to their trigono- 
metrical functions. On paper, owing to the algebraic sim- 
plicity of equation (5), the system appears to be less ccmplex 
than before; in fact, the simplification is illusory, for it 
actually involves at least three references to tables tor each 
surface, and these take up most of the time and also involve 
most chance of error. The system as given above is the form 
that is generally used ; but experience has shown that it is 
quite as quick to avoid reference to tables at all. The calceu- 
lations are generally arranged in three parallel columns as 
follews :— 


Central 
Angle. Sine. perpendicular. 
yy —> sing, 
Oe <o sin 9, ae 4 1 
1 oe sin 9, ars hy 
¥ 
Yi 7 5A sin py —— hy—h,! 
Y 
Pe ore SIN. Oe er hs 
1 
Oo <— SIN 05 <—. h,f 
Y 


The arrows show the order in which the operations are carried 
out. 
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_ It frequently happens that, after the ray has been traced, 
it is required to find the positions at which it is intersected by 
the focal lines formed by a pencil of neighbouring rays. The 
calculation of the radial focal line was considered in the previous 
Paper. It is most simply effected by finding the length of the 
ray intercepted between successive surfaces and the power 
of the various surfaces for that ray. The formula for the inter- 
cepted length d, corresponding to the axial length t is 


A, =A4,C08 Yr+7,COS 9’, —Tr+1 COS O41 iu. "(0 
and the power of the surface for this ray, Ky, is given by 
Ki, ==(UxCOs O— fy—1 COB), . «- , - (8) 


The position of the radial focal line is found by regarding d, and 
K) ast, and k,, and applying the formule for paraxial rays to 
find conjugate points. The determination of the transverse 
focal line is considered later. 

With the above treatment, as with all other systems that the 
author has met, special formule are necessary when flat or 
nearly flat surfaces are encountered. The application of the 
formule as they stand to nearly flat surfaces would result in 
considerable loss of accuracy, because two successive a’s are 
very great, but of opposite sign, and with a given number of 
. figures in sin py the sum of the two successive products a sin y is 
uncertain to an extent which increases with the magnitude of 
the radius. To avoid this uncertainty, and yet be able to 
treat all surfaces alike, whether they are flat, nearly flat or 
very considerably curved, an alternative method of calculation 
has been devised. 

Consider for a moment the conditions which such a system 
must satisfy. In the first place, all points of reference must be 
at a finite distance from the portions of the surfaces operative 
in producing refraction ; thus, no reference is possible to the 
centre of curvature or to the point of intersection of a ray with 
the axis, since either may be at infinity. Again, the radius of 
curvature may not be used, for this may become infinite ; on 
the other hand, its reciprocal, the curvature, may be employed, 
since it is always finite or zero. More generally, no lengths 

measured along the axis may be used if high accuracy is 
desired, because these are so variable im magnitude. Transverse 
distances, on the other hand, vary within small limits fixed 
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by the apertures of the various lenses, and their use will tend 
to give uniform reliability at all surfaces. Lastly, if the 
formule involve fractions, the denominators must in all cases 
be essentially constant in sign and large in magnitude. It 
will be shown that all these conditions can be satisfied. 

Take, first, a ray which lies in an axial plane. The figure 
shows the incident ray NP refracted at P as the emergent 
ray PM’. Let C be the centre of curvature of the refracting 
surface, and A the point at which the axis of the system meets 
the surface. Let M, M’ be the feet ot the perpendiculars to the 
incident and refracted rays from the centre of curvature C, 
and N, N’ the feet of those from the vertex A. The diagram 
. Shows the case in which the various quantities entering the 
formule are positive. Light travels from left to right, and 


the radius of curvature is positive for a surface presenting 
1ts convex side to the incident light. The inclinations y, y’ 
of the rays to the axis are positive when the rotation of a 
straight line from the direction of the axis through an acute 
angle in the counter-clockwise direction will bring it into paral- 
lelism with the ray. Similarly, the angles 9, 9’ of incidence and 
emergence are positive when a rotation through an acute angle 
in the counterclockwise direction from the direction of the 
normal to the surface at the point of refraction brings the 
line into coincidence with the direction of the ray. AT and 
AT” are perpendiculars from A to CM and CM’. The radius 
and the central perpendiculars will be denoted as before by 7, 
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hand h’, and the perpendiculars AN, AN’ by h and h’. The 
chord AP may be denoted by I. From the figure 


CM=r sin 9, 
CT =r sin y, 
AN=CM —CT, 

or h =r sin o—r siny. 


In this form the equation may not be used, because it involves - 
r. Dividing throughout by r and rearranging, the equation 
becomes 
simeo=sinwthR 1. | 2. « (9) 
a form which is free from objection. If, then, the incident ray 
1s given by yp and h, equation (9) enables the angle of incidence 
to be found. The angle of emergence is given by the law of 
refraction — 
Ww’ sin o’= Sin @ 
and the value of y’ by 
ol ea ek a 
At the next surface the length of the perpendicular from the 
vertex to the incident ray is evidently 
h’+¢ sin y’, 
where ¢ is the axial distance between the two vertices. It 
therefore only remains to determine the value of h’. This 
may not be found from the equation, 
h’=AN’=CM’ —CT’=r(sin 9’ —sin y’), 
since this involves r. Neither may r be eliminated by (9), 
pecause the denominator in 
sin 9’ —sin y” 


h’ =h-= 2 
sin @ —sin py 


becomes very small for long radii. Since, however, 9’ —y’ 
=o—y, this equation may be written 

h’_cos $(9’+y’)_cos 9’+c08 yp" 

h cosh(o+p) cos e+tcosyp’ 
and in this form the conditions laid down above are satisfied.* 
As the cosines of 9, 9’, y and y’ are required in dealing with 
the refraction of neighbouring rays, these are the quantities 


*It may be noted that h /cos }(+¥)= h ‘/cos 4(o/+y) =H, since 
the angles PAN, PAN’ are equal to $(9+¥) and 3/(9’+ ’) respectively. 
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that will normally be used. The equations of the new system 
are thus 


siti O\=sin ti ate ea a 
/8iN 05 = xsi SIO.) i 
Wr=O'r =O. Vned ar ao ne et (5) 

h Pmt Cosmos ae ee, eee 


COS Q,-+- COS Wy-1 


h NSS h +s Wr . . ° : 5 5 (12) 


The calculations are arranged very similarly to those set out 
diagrammatically above. Four columns are used—for the 
angles, their sines, their cosines and the perpendiculars from 
the vertices. It may be noted that the quantities found 
suffice to check all the calculations, since both the sines and the 
cosines of all the angles are known. In particular, absolute 
checks are available for all the quantities abstracted from the 
tables. 

The power of the surface / for the ray considered is derived 
as before from (8), but equation (7) has to be replaced. Several 
formule are available. If Td has been determined for the 
various surfaces, projection upon the axis gives 


d,cos wast, —$H2,Ry+hH 441 Ry41 ae he) odes (13) 
or projection upon the ray leads to 
d, =t,c0s Wy = V (Hd 2,—h 2 b/g = 2xea) : (14) 


When Hi is unknown, perhaps the simplest formula which 
does not involve the determination of any new quan- 
tities is 


SIN Q41-+8IN py 
. ic 
COS © 41+ COSY 


, Sin 9y-+sin 
A ny : 
COS @,-+COS Pr 


d, =t,cos Wr h + h Nee (15) 
Other interesting formule are obtained by projecting on to a 
line bisecting the angle between the ray and the axis. A 
simple formula of this class is 


d,—h = (Hay sin $0.41 — Ha sin $9',) sec $y. 


From this it “ollows that the optical length of the ray between 
the feet of perpendiculars to the incident and final emergent 
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paths from the two external vertices exceeds the axial path 
between these vertices by 

eX Hy sin 9’, {sec dy,_1 sec $¢—, —sec Fp, sec $¢’y} 


; Hi, a Sina 
cos dy,_1¢08 $ wa cos $¢ cos $0 


4 
M 


{cos $(9’,++-pa) 


: —C083(G.+Yr-1)} 
(h4=Ra)ma sin 94 _ 


cos sy,—1 cos dy, cos $¢ cos $9’y 


- 
M 


ais (ha— Ram sin 9 ea 
V i(1-+cos ya—1)(-+ cos y,)(1-+-eos 9)(1-+cos 9’n)} 


~ A number of other forms can be found for this path difference, 
but if is unnecessary to pursue the subject further here. 
Before leaving the consideration of a ray lying in an axial 
plane, it may be noted that the two series of formule using the 
centres of curvature and the vertices of the surfaces respec- 
tively as reference points may be used in conjunction with one 
another. Reference to the centres of curvature would be made 
for short radi, and to the vertices for long radii. All the 
formule remain as before, except that, on changing from one 
system to the other, instead of a, or é,, the axial distance 
between the two reference points is used. For example, if the 
centre of curvature is the reference point for surface 4 and the 
vertex for surface A+}, equation (1) will be replaced by 


h ati =A +h sin Wr» 


_ (16) 


where by =th —1>. 


_ If at surface 2+2 the centre of curvature is employed again, 
the transition equation will be 


Ays2= ‘agit, 41 SIN Yaz 
where Cpe =bhartly+e 


The corresponding changes in the formule for d are obvious, 
and need not be set down here. 

The extension of the formule to rays which do not always 
lie in the same plane may now be considered. As in the earlier 
Paper, let the ray incident in the medium sa, on the surface A 
be determined in direction by the direction cosines J,_;, My_1, 
N,_1 ; and the refracted ray in the medium jy, by Ly, My, Ny. 
The direction cosines of the normal to the surface at the point 
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where refraction takes place are Ix, m, mm. The previous 
notation involving © and y may be retained also, so that 
Ly, =CO08 yy ; 
and Lys + MWy1m + Ny -1 =C08 © 
Lyh-+- Mym+ Nam =cos 9’y 
In the two dimensional case the perpendicular from the vertex 


to the incident ray made angles 5 +ya-1and 5 + 9 with the axis 


and the normal to the surface. In the general case this will no 
longer be true. Let the direction cosines of the two perpen- 
diculars from the vertex to the incident and emergent rays be 
Cn, ik ON 
and ex» tks Or 
respectively. The co-ordirates of the feet of the perpendicu- 
lars are thus 
Qa Aha gaha 
and alias hr phx 
where the vertex of the Ath surface is taken as origin. The 
point of refraction is 
(1 —h), —1yM5 —TxMa, 

and the lengths intercepted on the rays between the perpen- 
diculars and this point are y/(Hi{—hyx) and +/(HIi—h 4). 
Therefore 

2 2, ™l—hb)—-aha  nmthha  nmtgrha 
es amie My Ma 

(17) 
Multiply the fractions, both numerator and denominator, by 
€&,fx,9, and add. Since 
Lat hMyat+gNy-1=9, 
it follows that 
Tr {0x—(Ahthim+gam)} = ha 

or Oly + fyma+ gym =ea —h rR or hee (18) 
the equation which leads to the value of the angle of incidence 
in the two dimensional case. Again multiply the fractions by 
1+), m,, m, and add as before ; this gives 


V (HHA — hix)= —(e+ealatAmat-gam) h/(cos pra eos @) 
=(h 7, —2e, h,)/(cos y,1+¢08 g) . . (19) 
by (18). Also multiplying by Z,_1, M@)_1, Ny; and adding 
V/ (HA hi)=n(cos y—1—C08 g).. . . (20) 
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The elimination of ./(H1; — hi) between (19) and (20) leads to 

Cos? 9, =cos? y_i+2eqh,R,—hiRy . .. . (21) 
as an equation for the determination of a. The value of 
9x will be derived from 

fy SIN 9, =Ly—1 SIN Or 
and K, from 
Ky =( jy COS ) -— Ma—1 COS) Rr. 

When @ has been found Hj is derived from (19). Equation 
(20) cannot, of course, be used for this purpose. It is next 


necessary to determine y,. This could be derived directly 
from the law of refraction 


Hy COS Wr —Ha—1 COS Pr—t 


Hy COS Ox —My.-1 COS QA= 


h 
Hs C08 Yr fat COB Yam 
1—F HER} 


but it is preferable to find ./(HIx—h’) directly, and thence 
obtain y, using a relation which may be written down by 
analogy with (20)—viz., 

COS Pa==COS Ort Ra — ha)... <m. 2 (22) 
The reason for this procedure is of course that, under the con- 
ditions which have been laid down, cos ya can be found from 
the value of W(Hx—h%); but the process may not be 


reversed. To obtain an expression for 1/(H;—h’?), write 
the equation of refraction in theform 


Jarr(COS Pr—COS OX) = a1" (COS Yx_1 —COS ©) 
—(1—A)ry(s4, COS 9 — 1 COS &) 


or pa (HR —h =a (HR—h})—-EBRK. . . (23) 
The equation : 
& hy=FHER, —cosp,V/(HX—h4). 2 2 |. (24) 


derived from a relation similar to (17), completes the determina- 
tion of the three quantities h X, exh x and y, from the given 
values of hx, ea hy and y,_1. 
The length on the ray in medium ya between the perpen- 
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diculars from the vertices of the surfaces A and J-+1 is evidently 
t, COS Wa, and therefore 


/ 
@xai Rapier Ripe mee. yeE1 —f,h, 


COS W) =- 
o ies Tix M), 
Ln / 
91 Nagi —Gr Do 
Ny : 
or ex Rast, RAS ESN ge ee 
and hr awi=hZ—-2ehxh-+h sin?y, . . (26) 


are the equations of transference to the new surface. 

To complete the calculations on the lines laid down in the 
previous Paper, formule are required for d,. Hither of the 
formule given earlier, 

dy cos Wy =h —4tH, RASH Ryo, ast ok 0 (13) 
or d,=t, cos pr — (HX — hR)-a/ (HE g1— Re *aga)s (14) 
may be used for this purpose. The latter is usually more 
convenient. 

It should be noted that in the preliminary calculations to 
which the above equations relate, neither the sines ot the angles 
involved nor the angles themselves need be found. The equa- 
tions relate entirely to cosines with the exception of 

{4 SIN O,=/y—1 SIN o, 
and this equation is in practice squared and written as a 
cosine equation, 


—COS? @y). 


The use of cosines instead of sines would be objectionable in a 
system in which reference to tables is involved, since a know- 
ledge of the cosine of small angles does not determine the angle 
itself to high accuracy. This does not apply to the system 
outlined above, where the object is not to know the final angles 
of emergence, but to determine the K’s, d’s and the cosines of 
the various angles of incidence and reiraction. The position 
of the emergent ray is obtained from a secondary calculation, 
involving these quantities, as shown in the spies Paper. 

As in the two dimensional ease, the formule just obtained 
may be combined with those contained i in the previous Paper. 
The equations of transierence are easily found and need no 
discussion. It will be noted that many alternative forms can 
be given to the equations for performing the calculations con- 
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sistently with the conditions laid down earlier. The system 
adopted is to be regarded as an illustration of a possible 
arrangement, and is not necessarily the best of its class. 

The operations may be summarised as follows: Given 
His xh, cosy 1,’ it is required to find h%,43, &41h adap 
cos y, and also dy, Ky, cos 9, COs 9}. 


COS? 9, =cos? wy, —( hR, —2e, hy)R, gives cos @ (21) 
ux(1 —cos? ox) = 2,4 (1 —cos? 9) » cos 9’ (10) 
Ky =(UCO8 OX — GCOS O) Ry » BK (8) 
PrvV/ FAR — hh Q) = py_1vV/ (HR - bh Kee Gx Daw (23) 
COS WY, =COS +R / (Hi Fy 5) 55 COS Wr (22) 
& hy =SHRR, —c08 pr (HX h 9) » &hx (24) 
O41 Axa =e hy —h(1 —cos? yr) s€xt1 x41(25) 
hah 2-26 hy 6+4(1 —cos? y) 5 A*%.41 (26) 


d, =t,cos wy, —V/ (HX — h a ae Gt Be = cas Be d, (14) 


It is evident on comparing these equations with those 
obtained in the earlier Paper that the conditions imposed have 
somewhat lengthened the calculations. The addition will, 
however, be found to be very slight. The appearance oi the 
equations gives a misleading impression of complexity, but this 
is simply due to the notation employed. Two square roots 
have to be extracted to determine cos 9 and cos 9’, the same 
number as in the earlier system. All else is multiplication, 
division or addition. The trigonometrical equations of Von 
Seidel, which in appearance are remarkably simple, will Le 
found to take up much more time. The latter involve at least 
20 references to tables for each refracting surface, and the 
calculations when completed afford less information than can 
be derived from calculations on the algebraic system. 

The earlier Paper already mentioned discussed the calcula- 
tion of radial focal lines. It may be of interest to add the 
corresponding arrangement which has been found most con- 
venient in dealing with the primary or transverse focal lines. 
As in the case of the radial focal line, the arrangement of the 
work resembles that followed for para xial rays. Put 


Hie=hecorg, cong. «1; . ~ (27) 
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and let Win, aoe be defined by the equations, 
Xr 


M1 eee ete ae 
ee OM» ee ee 
OM 1.a+1__ COS ® OMia_ rh ie, 
OMrx+1 COSQ OMr fa 
where 
poets Mir=Yv Wvo=2. 


The same value of ¥1,n is obtained by starting the calculations 
at the other end of the system, using the formule, 
COS 


_. Orn 
Arn cos o, tin He OM A SAS ROO. 


OMra—1n_€OSQ OMrn ey 
=F = aj? 
Gna: WO05.O. OM niin ie 


OM nn 
OM n 
The quantities defined by these equations have properties 
analogous to the corresponding paraxial functions. Rays 
arising from a point at a distance p, along the incident ray in 
the first medium measured in the direction in which the light 
travels from the point at which the first refraction takes 
place will be retracted through a point on the emergent ray 
at a distance, p,, from the point of refraction in the last surface 

measured in the same direction, where 


COS, Oia T = 
= =. as n . . . . 0) 
pomp eet, en Nan (30) 


where ==[, Mnn—Mn» Mnt+i n=, 


and 


mire cos Q. oM1 Mes ) 

ie: 1( On OMn . [Mins gre ¢. (81) 
provided the rays considered lie in an axial plane and lie close 
to the original ray throughout their passage through the optical 
system. Also if a small linear object lying in an axial plane 
and normal to the original ray meets the ray in the point 
distant , from the first surface, its image formed by neigh- 
bouring rays in an axial plane will be normal to the emergent 
ray, will meet this ray in the point determined by pn, and the 
linear magnification will be 9. 
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Many other formule can be contrived for these calculations, 
and it is easy to put down generalised formule. It will be 
found that none of these are so convenient as those given above 
for general use, though in special cases they may offer distinct 
advantages. As an illustration, take the formule, 


P=k, 
Oy 
P ,=cos? 9Pi,-1+P, TP.’ 
OP, A+1 2 OP, dy 
OPS wag WOP. jae 


with the symmetrical reverse formule, 


/ Cr. n 
Py n=cos? oN Prin tPy OP 
OP eis eit 2 OP Gs d,—1 
es 


It may be shown that 
Py n=M 1yn COS 1 COS Py" COS Go COS Po’ ... . COS Pn COS Ons 


where the continued product on the right consists of the 

cosines of all the angles of incidence and emergence. Also 
OP in __ cos On OM 1,n 

Dn, OP,, / ee Qn ° a AST ia 

Thus the terms involving the magnification and its reciprocal 
in the formule for conjugate distances would, with the P for- 
mule, require to be multiplied by the above continued product 
of cosines, and the labour of calculation in a complex system 
would be much increased. 

It follows from the formule for conjugate distances that 
Win is the power of the system for a pencil of primary rays 
lying close to the given ray, just as K,,, is the power for a 
pencil of secondary rays. In the application of the formule 
for the primary rays it must be remembered that 9 in general 
will differ from G@ for the same object and image points in a 
system corrected for astigmatism, because in the case of the 
former the dimensions of the object and image are not usually 
measured in planes normal to the rays. The appropriate 
valuefor 2 is to be found by orthogonally projecting the actual 
dimensions of the object and image on to planes normal to the 
incident and emergent rays. 
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